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Deterministic Model for Simulating the Predation 
of Toxorhynchites rutilus rutilus on Aedes aegypti 


By Dana A. Focks,! Jack A. Seawright,! and Donald W. Hall? 


ABS TRACT 


A deterministic computer model detailing the interaction of the con- 
tainer-breeding mosquito Aedes aegypti (L.) and the predatory larva of 
Toxorhynchites rutilus rutilus (Coq.) is presented. Simulation runs involving 
the release of Tx. r. rutilus adults indicate that releases resulting in one larva 
per container are sufficient to reduce Ae. aegypti adult density 75% in 20 
days. The slow development rate of predator immatures enables control to be 
maintained for several months. Simulations of predator release and the use 
of adulticides indicate that it is possible to obtain zero adult densities. Final- 
ly, the model demonstrates that the most important parameter in determining 
the degree of control established is the distribution of predator eggs. KEY- 
WORDS: Aedes aegypti (L.), computer simulations, insect control (bio- 
logical), insect control (chemical), mathematical models, mosquitoes, 
Toxorhynchites rutilus rutilus (Coq.). 


INTRODUCTION 


Historically, mosquito control has involved 
source reduction and insecticides. Predicting the 
outcome of these measures is simple and straight- 
forward. Currently, additional methods of mos- 
quito control are being studied, some of which 
involve the interaction of two species. Difficulty 
in understanding the dynamics of the interaction 
between two species makes predicting the out- 
come of new methods complicated. Attempting 
to optimize a strategy that utilizes two different 
methods in conjunction, e.g., predators and in- 
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secticides, is particularly difficult. Therefore, 
computer-simulation models are becoming in- 
creasingly popular as a tool in the development 
and evaluation of new control strategies. 

The purpose of this paper is to (1) present a 
deterministic computer model examining the in- 
teraction of Aedes aegypti (L.) and the predatory 
larva of Toxorhynchites rutilus rutilus (Coq.) 
and (2) simulate the population dynamics of Ae. 
aegypti under different control strategies involv- 
ing the release of Tx. r. rutilus adults and the use 
of adulticides. Pertinent to this discussion is an 
explanation by Mertz (1970) of methods used in 
life-history analysis. The utility and rationale of 
simulation in developing control strategies were 
reported by Conway (1970), Haile and Weidhaas 
(1976), and Weidhaas (1974). 

While many of the parameter values for Tx. 
r. rutilus are from laboratory experiments con- 
ducted under conditions obviously different from 
those expected in the field, there remains a very 
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real value in the model-building exercise. Char- 
acterizing the interaction of the two species 
highlights areas where additional information is 
needed. Also, simulating with the present model 
gives us insight into how best to approach a 
large-scale release experiment where prey control 
is to be attempted. Finally, the preliminary 
model gives us a framework within which to in- 
terpret the resulting experimental field data. The 
model is written in FORTRAN; the program 
listing is given in the appendix. 


MODEL DESCRIPTION 


The model described herein may be called a 
compartment model (Miller et al. 1973). Each 
stage of each species is represented by a number 
of storage registers within an array (an array is 
a group of registers). The registers represent 1- 
day age classifications of the various stages; a 
particular array represents the age distribution 
of a particular stage or instar. Except for the 
Tx. r. rutilus larval stage, the durations of various 
stages and instars are determined by the number 
of registers within a particular array. The model 
is made to cycle daily by replacing the contents 
of the next storage register with the contents of 
the previous register multiplied by the daily sur- 
vival for that particular stage and species; the 
output of the terminal registers of one array 
(stage) is the input to the next array (stage). 
Most values of daily survival for both species are 


fixed during any particular simulation run (table 
1); the daily survival for Ae. aegypti first and 
second instars is a density-dependent variable. 
The length of larval life of Tx. r. rutilus is also a 
variable, it being a function of the amount of 
prey available. Details of these variables ate 
presented later in the paper. 

The data used for the modeling of Ae. aegypti 
immatures were largely derived from Southwood 
et al. (1972) and for adults, from Sheppard et al, 
(1969). These life-table and ecological studie 
were conducted in Bangkok, Thailand, in hopes 
of correlating Ae. aegypti population dynamic 
with the known seasonal incidence of dengue 
hemorrhagic fever. Surprisingly absent from both 
studies was any information on fecundity and 
ovipositional patterns; estimates for these param. 
eters were derived from data on Ae. aegypti in 
northern coastal Florida (unpublished data), 
The data for Tx. r. rutilus came from laboratory, 
outdoor-cage, and field experiments conducted 
by Focks and Seawright (1977). 

According to Sheppard et al. (1969), the main 
Ae. aegypti breeding sites in Bangkok are the' 
earthen or ceramic water-storage jars (about 100-. 
to 200- liter capacity) found in association with 
all types of housing. These jars usually contain’ 
water throughout the year and are replenished 
with rainwater, tapwater, or riverwater. South- 
wood et al. (1972) state that the water-storage 
jars number about 150 per acre, and about 53% 
of these are positive for Ae. aegypti immatures 
at any particular time. Interestingly, both stud- 


Table 1.—Parameters and their values for model computations 


1 Ree figure 1 tee vate 


Symbol Description Value | 
tia ih ; Aedes aegypti F : i | 
SA Dailyeadult#survivalanc scene Coen ee Deter ee 88 | 
EF INZeDINCHAY Consagdedonuaococod Sh boone onic No. eggs/oviposition ...... 93 
}P) Successful’ pupation neh seen eee ee eee PCue scree 90 
SE IDEM KY CECA IIIAANAL dap poodanguseooogsoncoddoousacouabasene pete: 98 
S Daily survival for 3d and 4th instars and pupae ......... DCtiencra: 95 | 
C Density-dependent coefficient for calculating SIt .................. 0.01 | 
SI Density-dependent daily survival for Ist and 2d instars .. pct ...... 0-100 
AEDIST Containers positive for Aedes immatures ................ JX po dowe 53 | 
4 Toxorhy mchites rutilus indatS rs | 
SAT Day wart SUEVIVall, fer ase cee Ee ur eee eee pet...... 88 | 
FT Recundityanseecn eee eee No. eggs/female/d...... 1 | 
SIT Daily immature survival (larvae and pupae) ............ JOO csco ce 99 | 
DISTBN Containers positive for Toxorhynchites larvae ........... DCun Eee 50-90 


ies reported remarkably stable densities and dis- 
tributions of all stages from season to season. 
Adult densities of 1,100 per acre, or 7.3 adults 
per container, were reported. Table 2 presents 
the development times for immatures and the 
average numbers of, immatures per container. 
Southwood et al. (1972) note that 20% of the 
embryonated eggs hatch without the flooding 
stimulus described by Christophers (1960) and 
that approximately 50% of the remaining eggs 
hatch with each subsequent reflooding. Since 
they reported no data on the frequency of jar 
flooding, our computer model simulates no-rain 
and daily-rain situations. A final feature of the 
population dynamics of Ae. aegypti in Bangkok 
is the density-dependent survival during the first 
two larval instars. The stability of the mean 
number of immatures per container appears to 
stem from a paucity of larval food (table 2); the 
number of new third-instar larvae is a conse- 
quence of the available food and not the number 
of newly hatched first-instar larvae. An algo- 
rithm for producing the observed larval densities, 
which is largely independent of oviposition or 
frequency of hatch, is detailed in figure 1. 

The modeling of Tx. r. rutilus, with several 
exceptions, is similar to that for Ae. aegypti. 
Figure 2 reveals that (1) the three last immature 
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Table 2.—Development times, average numbers, 
and percentages of immature Aedes aegypti 
per water-storage jar positive for Aedes aegypti 
in Bangkok, Thailand’ 


Mean . Mean No. Percentage 
Stage development individuals of immatures 
time (days) per jar per jar 
BS Sees art tence 2 0) 
Larvae: 
Ist- and 
2d-instar ... 5.2 20.0 51 
3d-instar =... 3.2 95 24 
4th-instar 6.5 8.0 20 
PUpae eee weer DEP, 2.0 5 
Motaleyee se iL 7cil 39.5 100 


' Southwood et al. (1972). 
” Excluding eggs. 


stages are separated into three arrays, (2) eggs 
hatch independently of rainfall on the second 
day after oviposition, and (3) oviposition begins 
6 days after eclosion and occurs daily thereafter. 
Figure 2 also shows laboratory-reared adult Tx. 
r. rutilus being released at 6 days of age. 

The nature and mechanics of modeling the 
predation of Tx. r. rutilus on Ae. aegypti are pre- 
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Figure 1.—Flow chart and algorithm for Aedes aegypti computations. Letters within parentheses refer to daily 
survival values listed in table 1. 
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FIGURE 2.—Flow chart for Toxorhynchites rutilus rutilus computations. Letters within parentheses refer to daily 
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survival values listed in table 1. | 
| 
( 
Table 3.—Larval development times of Toxorhynchites rutilus rutilus in the 
laboratory when fed Aedes aegypti immatures 
ny Average Average consumption of | Average Average consumption of 
aie development Ae. aegypti development Ae. aegypti 
eae us time Ist and 2d instars time 3d and Ath instars and pupae 
Hae (days) Total Daily (days) fi Total: ? Daily ; | 
1 ie? 6.4 5.3 ED, 16.4 15.3 
2 2.6 6.6 2.5 2.6 16.6 12.5 
3 8.0 260 7.5 3.8 8.0 alt 
4 32.0 2600 18.8 8.0 73 9.1 
Total <P" 243 Been BGO UN ame a be 1516 yO. " 
‘The numbers reported consumed are Ist- and 2d-instar prey, since early-instar 
predators cannot consume late-instar prey. 
’ Estimated from 3 replicates. 
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sented below. Table 3 shows the development 


times of and the numbers of prey devoured by 
various instars of Tx. r. rutilus when offered 
either early (first- and second-instar) or late 
(third- and fourth-instar and pupae) Ae. aegypti 
immatures. Table 4 shows how long the various 
\Tx. r. rutilus instars can fast before death, and 
table 5 presents the duration of development for 
larval instars when fed Ae. aegypti larvae or 
tropical fish food. 

The situation to be modeled is as follows: 
(1) Third and fourth instars of Tx. r. rutilus 
‘develop faster on a diet of late instars and pupae 
than on a diet of early-instar prey. (2) All larvae 
develop at a rate proportional to the amount of 
prey available until the amount of prey exceeds 
that given in table 3. (3) All predator instars can 
survive without prey for a period of time, and if 
provided detritus, can develop at a very slow 
rate up to eclosion. (4) The daily immature sur- 
vival (SIT) is constant, irrespective of diet 
(unpublished data). 

In the model, first- and second-instar predators 
(array T12) are assumed to eat before third- and 
fourth-instar predators (T3 and T4, respectively ) 
on any particular day. This assumption simplifies 
the programing of predation and does not, con- 
sidering the relative amounts of prey consumed 
by each instar, introduce significant errors. The 
number of storage registers within arrays T12, 
T3, and T4 (16, 32, and 96, respectively) do not 
correspond directly to 1 day each. This is because 
the larvae are moved incrementally within an 
array as a function of prey size and numbers of 
prey available. If no prey is available, all the 
predators are incremented one register forward. 


Table 4.—Larval life of Toxorhynchites rutilus ru- 
tilus in the laboratory when deprived of food 


Instar in No. Age at 
which starvation individuals beginning Daysibetore 
began observed of test (days) death 
1 9 3 0 6.9+0.3 
wy, 10 4or5 8.4+0.5 
3 10 9 or 10 18.0+6.5 
24 8 13 '59.0+22.4 


! Temperature of 27° C. 

“Temperature of 27° C during first 11 days, then 
raised to 30° C. 

‘One individual fasted 88 days. 


If the numbers of prey equal or exceed the num- 
bers in table 3, the larvae are incremented 4, 8, 
and 12 registers forward in arrays T12, T3, and 
T4, respectively. Intermediate numbers of prey 
result in incrementing the predator arrays by an 
amount proportional to the available prey. The 
Ae. aegypti immature arrays in containers posi- 
tive for Tx. r. rutilus (XA and YA) are updated 
daily for the effects of predation. Since cannibal- 
ism does not occur among the predators at the 
densities considered here, even when there is no 
prey, it is not modeled. 

Following a release of Tx. r. rutilus adults, one 
would expect, on the basis of the distribution of 
Ae. aegypti larvae (53% positive) and an as- 
sumed predator distribution of 80% (unpublished 
data), to obtain four types of containers—42.4% 
positive for both species, 10.6% positive for Ae. 
aegypti only, 37.6% positive for Tx. r. rutilus 
only, and 9.4% containing neither species. Since 
it is assumed that no Tx. r. rutilus adults result 


Table 5.—Duration of development for instars of individually reared Toxorhyn- 


chites rutilus rutilus at 28°+1 
tropical fish food 


C when fed diet of Aedes aegypti larvae or 


Time increase, 


Mean development time (days)! when fed— 


ae fish-food fed 


Instar 
1 ED, 
2, 2.6 
3 3.8 
4 8.0 
Total+standard deviation 


Ae.aegypli 
(n=48) 


15.6+1.4 — 


Fish food. 


over prey fed 


(n=83) (eb 
8.6 720 
5.2 200 
14.4 380 
79.3 990 
~ 107.5+19.8 690 


' According to the arithmetic method of Southwood et al. (1972), each entry is the 


time required for half the population to molt to the next instar. 


Adult Release 
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FIGURE 3.—Predator and prey distribution in containers. Percentages in parentheses refer to frequency of each type 
of container when 53%, are positive for Aedes aegypti and 80% are positive for Toxorhynchites rutilus rutilus. 


from oviposition into containers devoid of prey 
and since one type of container produces neither 
species, these containers are not specifically mod- 
eled in the program but are accounted for with 
scaling factors (fig. 3). To represent the salient 
features of the interaction in the field and to 
represent the resulting Ae. aegypti adult density, 
two containers, both positive for Ae. aegypti, are 
included in the model. Only one of these is inter- 
faced with the predator subroutine. Ae. aegypti 
eclosion and Tx. r. rutilus oviposition are scaled 
by factors of DISTBN and AEDIST, respective- 
ly, to depict the frequency of the two containers 
in the field. Because of the scaling factors, the 
numbers of Ae. aegypti and Tx. r. rutilus adults 
are in terms of one container, and multiplying 
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Scaling factors 


Tx. r. 
rutilus 


adults 
DISTBN 


12% daily 
mortality 


Both 
species 
(42.4%) 


47% 


Ise” io 
rutilus 
only 
(37.6%) 


the output of adults by the number of container 
per area gives an absolute population estimate. 


RESULTS OF ‘SIMULATION 

Figures 4 and 5 represent model-generated Ael 
aegypti adult densities per container for no-rain 
and daily-rain-into-every-container situations, re- 
spectively. Figures 6 and 7 depict the correspond- 
ing total number of Ae. aegypti immatures pel 
container for the two rainfall situations. In each 
instance and in all subsequent simulation runs, 
the model program was initialized on day 1 with 
7.4 1-day-old adult Ae. aegypti. The resulting 
numbers within each stage were largely inde- 
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FIGURE 4.—Aedes aegypti adult density per container with no control measures applied (no rain). 
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FIGURE 5.—Aedes aegypti adult density per container with no control measures applied (daily rainfall into every 
container). 
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FIGURE 6.—Aedes aegypti immature density per container (no rain). 
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FIGURE 7.—Aedes aegypti immature density per container (daily rainfall into every container). 
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pendent of the number and stages used in ini- 
tializing the model. The graphs were smoothed 
by plotting 4-day moving averages. The cyclic 
nature of the daily-rain situation in figures 5 and 
7 dampened with time, 

A comparison can be made from table 6 be- 
tween the number of Ae. aegypti individuals per 
stage per container observed by Southwood et al. 
(1972) and the number of individuals per stage 
per container generated by the model when no 
control measures were applied. The total number 
of immatures per container generated by the 
model exceeded the total reported by Southwood 
et al. This disparity was allowed to remain be- 
cause the difference was not large, especially in 


Table 6.—Number of Aedes aegypti per stage per 
container as reported by Southwood et al. 
(1972) and the number generated by the model 


Southw 


Stage 
ex et al. 


Model 
No rain Rain! Average 


ood 


Larvae: 
Ist- and 2d-instar . 20.0 
CYleinGiighe! coacoess 9.5 
Ath-inStareesce 8.0 
Pupae aa so00 PD) 
Motally Rerus rs 39.5 
PERT GUNS aN te een a meee FD 


30.0 70.0 50.0 
8.0 5.0 6.5 
8.0 5.0 6.5 
Dp 1.5 1.9 

482 81.5 64.9 
8.0 5.0 6.5 


1 Situation modeled is tha 
rain. 


t all containers receive daily 


Table 7.-Mean number of Aedes aegypti per container after contact adulticide 
application and after application when adult density exceeded seven or four 


adults per container! 


St Not Single 7 Adults 4 Adults 
abe OURGEIERO RE application” per container’ per container! 
ANiIt? tecrold ar eotaeetomey poe 8 8 3 D 
JinsrenehNKs soceoasppeaae 49 49 28 21 


1 See figures 8-13. 
“ Daily immature survival (SI) of 75%. 


* Daily immature survival of 84%. 15 applications in 262 days. 
‘ Daily immature survival of 88%. 21 applications during period 98-360 days. 
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Figure 8.—Aedes aegypti adult density following contact adulticide application 
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95% adult mortality with no residual effects. 
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“IGURE 9.—Aedes aegypti immature density following contact adulticide application on day 90 (no rain). 
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FIGURE 10.—Aedes aegypti adult density following contact adulticide application when Aedes aegypti adult density 
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the no-rain situation, and would have made the The effects on Ae. aegypti population dynam- 


model more conservative in subsequent control ics of contact adulticide application, which 
simulations. Density-dependent first- and_sec- caused 95% mortality among the adults with no 
-ond-instar daily survival (SI) for the no-rain and residual action, are presented in table 7 and fig- 
daily-rain situations averaged about 75% and ures 8—13. Figure 8 shows the overshoot of adults 
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FicurE 11.—Aedes aegypti immature density following contact adulticide application when Aedes aegypti adult 
density exceeded seven adults per container (no rain). 
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FIGURE 12.—Aedes aegypti adult density following contact adulticide application when Aedes aegypti adult density 
exceeded four adults per container (no rain). Spraying at this density required 21 applications in 262 days. 
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FIGURE 13.—Aedes aegypti immature density following contact adulticide application when Aedes ae 


density exceeded four adults per container (no rain). 
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FIGURE 14.—Aedes aegypti adult density following an adult predator release on day 98 that resul 


12 


application. Figures 10 and 12 show the increased 
frequency of application required to maintain an 
increasingly lower adult density. Because of this 
type of situation, population suppression has 
usually involved the use of larvicides in addition 
to adulticides (Gould, et al. 1970, 1971). 

Figure 14 shows the effects of a predator re- 
lease on day 98 that resulted in one predator egg 
in 80% of the containers during a period of no 
rainfall. A 75% reduction in Ae. aegypti adult 
density occurred within 20 days. The Ae. aegypti 
adult blip beginning about day 170 resulted from 
the eclosion of the predator. The effects of the 
subsequent second generation of predators re- 
sulting trom the release are shown after day 185 
in figure 14. 

Until actual field trials are made, the reality of 
the effects of a second generation will remain un- 
certain. The model does depict, however, smaller 
and smaller subsequent generations of predator. 
These results are consistent with known equilib- 
rium densities of predator and prey in nature. 

Figure 15 reveals that one predatory larva was 
sufficient to deplete a water-storage container in 
6 to 8 days. The model further shows that con- 
tinued oviposition by declining numbers of Ae. 
aegypti provided food, allowing the predator to 
develop and eclose in approximately 53 days. 

Simulations involving larger adult releases that 


produced greater numbers of predator eggs per 
container did not appreciably alter the Ae. 
aegypti decline. This is understandable because 
one predator per container was sufficient to stop 
prey breeding in that particular container; there- 
fore, the additional predators were superfluous. 
Larger numbers of predators per container fur- 
ther exacerbated the prey shortage, producing 
even longer development times. In the field, a 
variable number of predators per container would 
likely result in asynchronous predator eclosion. 

Since rainfall increased the number of young 
prey larvae, five predators per container were re- 
quired to eliminate Ae. aegypti breeding in con- 
tainers positive for Tx. r. rutilus. Figure 16 shows 
the failure of three predators per container to 
establish control. Furthermore, the increased 
food supply caused shorter predator development 
times, which resulted in shorter periods of control 
from the first cohort of predators. The ramifica- 
tions of this would depend on the results of sub- 
sequent generations of predators resulting from 
the original release. At any rate, the lower num- 
bers of predators required in the dry season to 
initiate control indicate a logical time to attempt 
control. 

Figures 17, 18, and 19 represent expected Ae. 
aegypti adult densities when the percentage of 


(Continued on page 16.) 
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Figure 15.—Aedes aegypti immature density following release of one predatory larva per container on day 98 in 80% 
of the containers positive for Aedes aegypti (no rain). The blip at 160 days is a result of predator eclosion from 


the container. 
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FIGURE 16.—Aedes aegypti adult density following an adult predator release on day 98 that resulted in three preda-| 
tory larvae per container in 80% of the containers positive for Aedes aegypti ( 


FIGURE 17.—Aedes aegypti adult density following an adult predator release on day 107 that resulted in one predator) 
egg per container in 60% of the containers positive for Aedes aegypti immatures (no rain). 
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FIGURE 18.—Aedes aegypti adult density following an adult predator release on day 107 that resulted in one predator 


egg per container in 80% of the containers positive for Aedes aegypti immatures (no rain). 
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FicurE 19.—Aedes aegypti adult density following an adult predator release on day 107 that resulted in one predator 


egg per container in 90% of the containers positive for Aedes aegypti (no rain). 
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containers receiving one predator egg per con- is the distribution of predator eggs (—DISTBN). 


tainer is 60%, 80%, and 90%, respectively. The It seems likely that energy would be well spent 
preceding discussion and these figures lend sup- in improving predator distribution through re- 
port to the idea that the most important param- lease of smaller numbers of adults at more sites, 
eter determining the degree of control established Simulations involving the use of adulticides 
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FIGuRE 20.—Aedes aegypti adult density following contact adulticide application on day 90 and subsequent predator 
release on day 98 (no rain). Compare results with those in figures 8 and 14. 
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FIGURE 21.—Aedes aegypti adult density following a predator release and five subsequent contact adulticide applica- 
tions on days 98, 101, 109, 113, and 117 when Aedes aegypti adult density exceeded four adults per container 
(no rain). Compare results with those in figures 12 and 14. 
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and predators revealed the following: (1) One 
adulticide application shortly before or after a 
predator release did not significantly decrease 
the resulting prey density nor increase the rate 
at which it was achieved (fig. 20). (2) Several 
applications 2 or 3 days apart after a predator 
release resulted in immediate and sustained prey 
densities near zero (fig. 21). These results reen- 
force our assumption that maintenance of control 
is a function of the efficacy of the resulting 
subsequent predator generations. 


CONCLUSIONS 


As stated previously, the computer model is a 
useful tool in evaluating and developing new con- 
trol strategies. The utility of Tx. r. rutilus in 
controlling container-breeding mosquitoes will 
have to be demonstrated in the field rather than 
in a computer. Demonstrated here is the feasi- 
bility of control and a justification for larger scale 
fieldwork than has already been done. 

Field data on predator distribution, the effects 
of second- and third-generation predators result- 
ing from an initial release, and _ field-survival 
values for the immature predator stages will 
make it possible to develop an accurate stochastic 
model. This model, in turn, will facilitate the 
evaluation of a strategy involving insecticides, 
predators, and the sterile-male technique. 
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APPENDIX.—FORTRAN LISTING 


REAL GRAPH( 3000), DAY(360)> So ee oe bes SUMTA( 
1SUM1 (360) + SUM2(360)+ SUM3(3€0)>. piseons KX 2 ( 
1X30 260)-5 XK4( 360), X5(3€0) > X6( 360) « X7(360) 

DATA DAYs SUMAs SUMTAs SUM1L+ SUM25 SUM35 X1ls X25 
1X35 X4s X5Ss X6s X74 360*00e 360*00r 3600*00% 360*0 os 
1360*0+s 360*050s 36O0* 00s 360* 002 360* 050s 360*00es 36O0*00s 369*0e» 
VZECO*xO. “Ue , 


e) 
Ov 
° 


B2(3€0)s B3(360)s 8B4( 360) 


B1¢ Vs 
DATA B1+.H2.B32B4/ 360*0 « » 360 *00s 360* Oe » 300* 00/7 
REAL C1 (360) +C2(360)+C3(360) »C4(360) 
DATA ClsC2sC3sC4/7 360* Ce es 3EC¥KCe »360*00 s 360% 00/7 
REAL 0O1(360 ) 
DATA D1/360*0 0/7 
REAL 85( 360). C5( 360) 
DATA 8S. cS 7260 *002=60*00e/ 
RFAL A(B)sTAC7) 

A € TA REPRESENT AEDES & TOXe ADULTS 


REAL E(6)+X(5)eY(12) 
AEDES IMMATURES IN CONTAINERS WITHOUT TOXe 
WHERE E=EGGSs X=STAGES 162%6 Y=STAGES 3:s4EPUPA 6 


REAL TEMP12(24)5 TEMP3(40).5 TEMP4(108) 

DATA TEMP12s TEMP3s TEMP4 /24%* 0ee 40*0 05 108*0O-e 7 

REAL FT(6)%XA(5) +YAC 12) 
AEDES IMMATURES IN CONTAINERS POSITIVE FOR TOX. 
ET=AEe«EGGSe XA=AEeCSTAGES 162 & YA REPRESENTS STAGES 3-PUPA 


RE AL T12(°024)+.T3(40),T4(108)5TP(8)5 TEC 3) 
TOXORHYNCHITES IMMATURE ARAYS; T12=STAGES 162s T3=STAGE 3 
T4=STAGE-4-e- TP=PUPAs ETE=EGGS 


DATA ArsTAS B¥*0ee7*0e 4 

CATA EsXeV¥F 6*0025%*¥0e212*00 / 

DATA ETsKAsVYAS 600 25*00912*06 / 

DATA TLEesT3ZsT4eTPSY 24% 00 s40* 00 108 X00 rs B*Oe S 


DATA TEY 3*06 i, 
=0 
AEDIST=PROPORTION OF CONTAINERS POSITIVE FOR AEDES LARVAE 
SA=0.-38 
SA=ADULT DAILY SURVIVAL 
r= oS 
F=FECUNDI TY (EGGS/BATCH) 
P=0.90 
P=PROPORTION GF PUPATION WHICH IS SUCCESSFUL 
SE=0.e98 
SE=DAILY EGG SURVIVAL 
S=0.e95 
S=DAILY SURVIVAL FCR STAGES 394 & PUPA 
C=0.01 


G€=DENSTITVY DEPENDENT COEFFICIENT FOR CALCULATING 
SI»WHERE SI=DAILY LARVAL SURVIVAL sSTAGES 1 E26 


***XTOXe PARAMETERS ¥ Xx 
SAT=0-¢ 88 
SAT=DAILY ADULT SURVIVAL. 
i160) 
FT=DAILY EGG FRODUCTIONe 
SIT=0.99 
SIT=DAILY IMMATURE SURVIVAL» STAGES 1-—-PUPA. 
OILSTBN=PROPURTION OF CONTAINERS POSITIVE FCR TOXe EGGS. 
CISTBN=0.80 
eK kok kok ak ok kk kk gk ok eck ke ak ak OR ok OK kg ok Kk kok kk ak dk ak kk ke aK OK i ake ak de a ko ok ok Ac 
ITRATN=1 
I=I1I+1 


eee iar AECES wITkF 764 ADULTS/SCONTAINER 
TeFQel) AC 1)=7 
-CW REPRESENTS AN ADULTICICE APPLICATION CAUSING 95% MORTALITY 
TetTe-98) GOTO 1501 
DUET DENSITY WEEN TUEV TS APPELED 


19 


1501 


‘evane) 


150 


20 
30 


ane) 


100 


200 


(e@QH@) 


151 


20 


R=1000 
SA=O¢ BB 
IF (SUMACT=1)e¢GTeR) SA=0205 
TF (SAeEGQee05) GOTO 1500 
GOTO 1501 
WRITE (621502) I 
FCRMAT (1XeI13) 
CONTINUE 
IF (LeEGe1l07) TE(1)=1 
**KKAF DES ADULTS**x* 


A(8)=A(7)*SA 

A(7)=A(6)*SA 

ACE) =A(5)*SA 

A(S)=AC4)*SA 

A(4)=A(3) ¥*SAtA(8)*SA 

A(2)=A(2)*SA 

A(2)=A(1)*SA 

A(1)=(Y(12) *¥P*¥(1 -DISTBN))I4(YAC12)*P*DOITSTBN) 
***KAEDES LARVAE IN CONTAINERS WITHOUT TOX>s’ LARVAE* xx 

FOLLCWING MOVES STAGES 3-PUFA 

YC 12)=YC1T) *S 

Y(11)=Y(10)*S 

Y(1C)=Y(9)*S 

Y(S)=Y(8)*S 

Y(8)=Y(7)*S 

YO 7I=AVYCEV*S 

Y¥(6)=Y(5)*S 

Y(5)=Y(4)*S 

Y(4)=YC3)*S 

Y(2)=Y(2)*S 

Y(2)=Y(1)*S 

FOLLOWING CALCULATES SI FOR STAGES 1 AND 2 

SUM=X (1)4X02)4+K°03) 4X04 4X05) 

IF (SUM-EGe0) SI=1020 

IF( SUMefQe0) GOTO 30 

TF (SUMeLT 215) GOTO 10 

IF (SUMeGTe120) GOTC 20 

XNT= SUM*(EXP(-C*SUM) ) 

SI=XNT/SSUM 

GCTC 30 

S1=0.88 

GOTO 30 


NG MOVES STAGES 1 AND 2 
y*ST 


ee ee en 
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Ll 
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wwveuewf 
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HE FOLLOWING CONTRCLS EGG HATCHE OVIPCSIT ION 
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AEDES OVIPOSITION INTO TOX-FREE CONTAINERS 
E(1)=A(4)*F*005 

EGG HATCH AS FUNCTICN OF RAINFALL (TRAIN) 
IF (IRAINeEQel1) GOTO 100 
See a 


Pd 
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au 

6) *O00e5+(E(5)*00e ) 
ee i 
5) 
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AEDES LARVAE IN CONTAINERS WITH TOXe LARVAE *** 


FOLLCWING MOVES STAGES 3-PUPA 

VAC LEI=AVYACLII*S 

YACLLI=VYAC1O)*S 

YAC10 

YA(9) 
) 


emt tl Gr el Ga) PS SP IK SR} 


11 
21 
Sys 


101 


201 
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al 


ol 


29 


De 
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41 
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HNWFUO 
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CULATES SIX FOR STAGES 1 AND 2 
SCMA=XA(1) (2)#XA(3)+4+XA(4)4XA(5) 
IF (SOMAcEQ SIX=120 
IF (SOMA2EQe0) GOTO 31 
IF (SOMAetTe15) GOTEG 11 
IF (SOMAeGTe120) GOTO 21 
XNTA= SOMA* (EXP(-C *SOMA )) 
SI X=XNTA/SOMA 
GOM@) Sil 
SIX=0.-88 
GG To) 31 
SIX=120/S0MA 
FOLLOWING MOVES STAGES 1 AND 2 
YAC1)=XA(S) *SIX 
KA(5)=XA(4) *SIX 
XA(4)=XAC3)*SIX 
XAC3)=XAT2)*SIX 
KA(C2)=XAC1)*SIX 
THE FOLLOWING CONTROLS EGG HATCHE OVIPOSITION 


Ox PDH He eH HH 
~rpFnnnnunn 


BiIt6) — CE (5)+E 7 (6) *Se& 
Ex SE C4 SE 
ET(4)=ET( 3) *SE 
EeyCsySseu (2). 
EyGe)=EuCl) *Se 


AEDES OVIPOSITION INTO CONTAINERS POSITIVE FOR TOXe 
ETCLI=A(4)*F* 005 

EGG HATCH AS FUNCTION OF RAINFALL (TRAIN) 
IFCIRATINeEQe1) GOTO 101 
XACLI=ETCS)*O062 
ET(S)=ET(5) #028 


GENE, 207 
KACLI=ET(6)*O054+(ET(5)¥* 002) 
ET(6)=ET(6)¥*005S 
ET(SI=ET(S) *008 
CONTINUE 
*x*AEBELOW REPRESENTS TOXe IMMATURES € THEIR PREDATIUN*** 
ON AEDES LARVAE IN THE ABOVE SECTION 

THE FOLLWCING INTERFACES 1&2-STAGE AEDES 

WITH 1€2-STAGE TOXe 
SMXA1=0 


SMXA1=NO5co 1E2-STAGE AEs LARVAE IN AEo/TXe CONTAINER 
DS 13 KN=1+5 

SMX A1=XACKN)+SMXAI1 

SMT12=0 

SMT12=NGe 162-STAGE TOx*e LARVAE 
DC 71 M=14.16 

SMT12=SMT12+T12(M) 

IF (SMXA1e¢GTeSMT12*4) GCTO 29 

ITF (SMKA1*e¢EQeO) I[Z=1 

IF (SMXA1¢EQe0) GOTO 39 

IF (SMT12e¢EQe0) GOTO 601 

IZ= INT(SMXAIL/SMT12) 

Mle <o1G/e Oe li Z—i 

CC 61 L=145 


XAC(L)=0 
GOTO 39 
IZ=4 


ATE=(4¥*SMT12) 
NECA its i= le10))) (GOTO) sg 
BE=ATE/SMXA1 

CCc=1-88 

DERS2) K=145 
KACK)=XA(K)*CC 
CONTINLE 

DC 41 J=1516 
TEMP12( 441 Z)=T12¢05)*SIT 
CONTINUE 

OC 999 MM=1524 


ANNDOA 


fa) 


ANN 


(A) 


999 T12(MM )=TEMP12(MM) 
601 CONTINUE 


THE FOLLOWING INTERFACES 3€4-STAGE TOX. WITH ALL STAGES 
OF AEDES IMMATURES»e MGDEL ASSUMES AEDES 1E2-STAGE LARVAE 
ARE EATEN FIRST FOLLCWED HY 3,46P-STAGE AEDESe NUTICE 
THAT FOR EACH DAYsTHE 1€&2-STAGE TOXe EAT FIRST. 
T301)=712(17)4T12(18)4T12(19)4#T12(20) 


DG 400 K30=17»20 
T12(K30)=0 
400 CONTINUE 
SMXA2=NO0 AEDES LARVAE sSTAGE 1LE2eAVAILABLE AS FOUD e 
SMXA2=0 
DC 53 N=15 
53 SMXA2=SMXA24XACN) 


SMY=NO56 3546 PUPAL-STAGE AEDES AVAILABLE AS FOOD. 
SMY=0 
DO 22 LL=1.12 
32 SMY=SMY4YA(CLL ) 
TOTU=SMXA24 (SMY*7e1111 ) 
NOTE-1 334 OR PUPA-STAGE AEDES IS EQUIVLENT TO 7.111 
1 OR 2-STAGE AEDES LARVAE.’ TOTU=AMOUNT OF FOQD AVAILABLE. 
ST3=NO5o 3-STAGE TUXes ST4=NO0 4-STAGE TOX. | 
ST3=0 
DC 42 I1=1,32 
42 ST3=ST3+T3(11) 
ST4=0 
DO 85 JJ=1296 
ST4=ST44T4( JU) 
85 CONTINUE 
REQ3E4 ARE THE NOe CF AEDES EQUIVe NECESSARY TO MOVE 1 
REGULAR DAYe 
REQ3=ST3*14 
REQ4=ST4*64 
TREGQ=REQ34REQ4 
IF (TOTUCEQsO) IT3=1 
IF (TOTUcEQe0) IT4=1 
IF (TOTUCEGe0O) GOTO 105 
IF (TREQeEQe 0) GOTO 108 
IF(TRFEFQeLEeTOTU) GOTO 115 
DC 22 KKk=1 55 | 
XA( KK )=0 | 
22 CONTINUE | 
DG 12 Vera, he | 
YA(KI)=0 | 
12 CONTINUE 
Q=TOTU/TREQ 
IT3€4 INCREMENT THEIR RESPECTIVE ARAYS WHEN LESS THAN MINe FOOD 
IS AVAILABLE. 
IT3=INT(@*Q) 
IT4=INT(12*Q) 
IF (IT3eEQe0) IT3= 
IF (IT4eEQe0) IT4= 
GOTG 105 | 
115 IT3=8 
IT4=12 
ITERATIONS FOR 3&€4-STAGE TOXe WHEN AMPLE FOOD AVAILABLE. 
BAL=TREG-SMXA2 
BAL HERE REPRESENTS THE NOe UNITS REQDe FROM YA, 
IF (BAL eGTeO0) GOTO 90 
SIAX=ABS(BAL/SMXA2) 
DC 91 MI=155 
91 xXA(MI)=XA(MI)*STAX 
GCTO 105 
90 CONTINUE 
DO 92 MJ=195 
92 XA(MJ)=0 | 
BAL=ABS(BAL) 
SIYA=((SMY*7e111)-BAL) /(SMY*7e111) 
DOL 2'S) SUN=Wenise | 
125 YACIN) =YACIN)*SIYA 
105 CONTINUE 
DC 106 KJ=1,32 } 
TEMP3(KJFtIT3I=T3(KI)*SIT 
106 CCNT INUE 
DC 998 M5=1,40 


T2(M5) =TEMP3(M5) 
$98 CCNTINUE 

DC 107 KM=1.96 

TEMP4(KMFtI1T4)=T40KM)*SIT 
107 CONTINUE 

DO 997 M6=1.108 

T4 (M6 )=TEMP4(M6) 
997 CONT INUE 

T4(1)=0 

0G 555 LM=33+40 

T4(01)=T401)4+T3(LM) 
555 CONTINUE 

NO 401 K31=3340 

PSICKSi19) —0 


S396 CONTINUE 


402 CONTINUE 


TA(7IY=(TACT7 )#TAC6) )¥SAT 
TA(6)=TACS)*SAT 
TACSI=TAC 4) *¥SAT 
TA(4)=TA(3)*SAT 
TA( 3)=TAC2]) *®SAT 
TA(2)=TAC1)*SAT 
WAG I — MP7) * Sit th 
TE( 2)=TE( 2) 
ENGen) — sels) 
TOXe OVIPOSITION & EGGS 
TE(1)=TAC7) *0.5*F T¥AEDIST 
Teli2 ee) —neES)) 
GUC P= Vi Gl) VAC 2 RY. CSCS) 
C2CTI=AVCS)4+V(6)4V07)4+Y(08)4Y°9G )tY( 10) 
SSC WSN Weep) 
Ciec &€3 ARE THE NOe OF AEe STAGE 3934 AND PUPA 
IN CGNTAINERS WITFOUT TXe 
C4¢(1)=SI 
cS5(1)=S1ix 
SUMA & SUMTA ARE THE NCe OF AEe AND TOXes ADULTS. 
SUMA(I)=0 
(oye) fsioh eht=sat oie} 
SUMACI)=SUMACT) +A J1) 


801 CGNT INUE 
SUMTACT)=0 
DONS OZ I2— Went 
B02 SUMTACTI )=SUMTACT)+4 TACJ2) 
DAY(I)=I 
D1(T)=0 


803 D1(I)=D1 
X01) =0 


204 X2(1)=Xe 


I) 

MMATURFES IN CONTAINERS WITHOUT TXe LARVAL 
X3(1I)=0 
OC 805 J6=135 

205 X3(T)=X3(1)+XA(J6) 
X4( 1 )=0 
DC 806 J7=15,12 

F06 X4(1)I=X4(1)+4+VA(JI7) 
SUM2( 1 )=X3(1)+XK40I 
SUM2=SUM OF AEs IM 
xXS(I)=0 
DE 807 J7=12916 

B07 XS(I)=xXS(1)+h1l2( 37) 


) 
MATURES IN CONTAINERS POSITIVE FOP TXe LARVAE o 


fa) 


24 


8028 


809 


308 
304 
305 
306 
307 
501 
502 
503 
504 
505 


471 


700 


701 


702 


705 


300 


XE(T)=0 

DC 808 J8=1,32 

xE(1)=xX6(1)+T3(48) 

X7{ 1 )=0 

DC 809 J9=1.,96 

X7CTI=X70C1)4T4(0599) 

SUM3(0 1 )=X5S( 1) tX60T D+X7¢ 1) 

SUM2=SUM OF TXe LARVAE 

IF (let T2360) GOTO 1000 

DG 308 [=1,.,357 
X1(C7T)=(001¢01)4010141)401¢142)+01(014+3))74 

DO 3204 [=13,357 
BLCIy=CC1CT)+C101+1)4C10 [+2 )4+C101+3)) 74 

DO 305 [=1,357 

B2C1T)=(C2c 1) 4+C2C +1) 4+Ce(l+2)4+Ce2(1+3))74 

OO 306 [=1,357 
B3ACTI=ACC3(TI+C3CT+1)4C3(142)4+C3(14+3))74 

DO 307 I=1,.,357 

BOTS CeO GEOL EY RCA CT 2) ECAC 3) 74 

OC 501 [=1,.357 

SUMA( 1 )=(SUMACI4#1)#SUMA( 1) #SUMACL+2)+SUMA(143)) 74 
DQ 502 [=1,357 
SUMTACT)I=CSUMTACT)+#SUMTA( I +1 )4+SUMTACT#2)4+SUMTACI+3) 
0C 503 [1=1,357 

SUM1 (CT )=CSUM1 CI )+SUM1C0T41) +SUM10CT+2)+SUM1(0143))74 
DO S04 [=1,357 

SUM2 (I )=(SUM2( 1) 4+SUM2C I +1)4+SUM2(1+2 )+SUM2(14+3))74 
DC 505 1=1.357 

SUM3(1T )=(SUM3 CI )4SUM3 C141) +SUM3( 142 )+SUM3(14+3))74 
COC 471 1[=1.357 

PSCII=HCCSCTI+ACHOCT +1) ¢C5( 1+ 2)+C5(C1+3))74 

EFELOW PLOTS AEe ADULTS 

CALL PLCT1 (3951 0100111210) 

CALL PLOT2 (GRAPH. 36006 09 Vexs20en Ve) 

CALL PLCT3 (1H*»DAY(1 ) » SUMA(1),5,357) 

WRITE (€s700) 

FCRMAT (C1H1s "AEDES ADLLTS® ) 

CALL PLCT4 (12s12HAEDES ADULTS) 

WRITE (651200) 

BELOW PLOTS TXe ADULTS 

CALL PLOT] (3251 0100111510) 

CALL PLOCT2 (GRAPH. 3600 9002200 20e) 

CALL PLOT3¢(1H*®,DAY(1)2SUMTA(1),357) 

WRITE( E0701) 

FURPMAT (1H1»s *TOXORHYNCHITES ACLLTS®) 

CALL PLCT4(11,11HTOXe ADULTS) 

WRITE (691200) 

BFLOW PLOTS AEe IMMATURES IN CONTAINERS WITHOUT TXe 
CALL PLET1 (3451210011110) 

CALL PLEOT2Z( GRAPH » 2E€00 29 Co 014005 00 ) 

CALL PLOT3¢(1H*,DAY (1 )o SUMI1(1)5357) 

WRITE (65702) 


FORMAT (LHe *AEe ITMMATURES IN CONTAINERS WITHOUT TXe 


CALL PLCT4(155915 HAEDES IMMATURES) 
wkIT& (651200) 

BELOW FLOTS AEe STAGE 1 & 2 
CALL PLCTL (3551210011110) 
CALL PLECT2 (GRAPH >. 32606 29 002 140 
CALL PLCT3 (1IH*,DAY( 1)» X1( 1). 
WRITS (€5705) 


FORMAT (1H1s *AEs STAGES-1E2 IN CONTAINERS WITHOUT TXe LARVAE ® ) 
CALL PLCT4 (9,9HAECES 182) 

WRITE (€01200) 

BELOw PLOTS AEe STAGE = 

CALL PLCT1 (3251 0100111510) 

CALL PLET2 (GRAPH s36060 2002500200) 

CALL PLOTS CLIH*.DAY(1)2B1(1)%2357) 

WRITE (65300) 

FORMAT (1H1s*AEe STAGE-2 IN COCNTAINERS WITHOUT TXe LARVAE ® ) 
CALL PLCT4 (13513HAEDES STAGE 3) 

WRITE (621200) 

BFLOW PLOTS AEe STAGE- 4 

CALL PLCT1 ¢€3,510100e111510) 

CALL PLOCT2 (GRAPH .»3606¢« 200x250 ) 

CALL PLCT3 (1H*.DAY( 1) »B2(01)+2257) 

WRITE (6.301 ) 


44 


IMMATURES® ) 


301 


302 


303 
x 


450 
1 


703 


704 


1200 


SENTRY 


FORMAT (1H1Le*AEs STAGE-4 IN CCNTAINERS WITHOUT TXe LARVAE ® ) 
CALL PLOT4 (13s13HAEDES STAGE 4) 

WRITE (691200) 

BELOW PLOTS AEe PUPAE 

CALL PLOT! (3651%102111510) 
CALL PLCT2 (GRAPHs3606¢ s00e2106 
CALL PLCOT3 (C1H*esDAY(1).B3(1)> 


WRITE (62302) 

FORMAT (1H1e*AEe PUPAE IN CONTAINERS WITHOUT YTXe LARVAE * ) 
CALL PLOT4 (1llellHAEDES PUPAE) 

WRITE (631200) 

BELOW PLOTS SI 

CALL PLOTL (3251 21002111410) 

CALL PLGOT2 (GRAPHs 36064 200 sles0e) 

CALL PLCTZ3 (1H¥*¥»sDAY(1)284(1)2357) 

WRITE (65303) 

FORMAT (1H1lse* IMMATURE SURVIVAL (SI) FOR STAGES 16€2 AEDES iN 
CONTAINERS WITHOUT TXe!*) 

CALL PLOT4 (2s2HSI) 

WRITE (631200) 

BELOw FLOTS SIX 

CALL PLGT!1 (€329512100111.10) 

CALL PLOT2 (GRAPH, Se ag aa Oe) 

CALL PLOTS (1H*.» DAY(1) >» SHO issow) 


WRITF (62450) 

FORMAT (1H1»es* IMMATURE SURVIVAL (SIX) FOR AEDES STAGES-18&2 IN 
CONTAINERS POSITIVE FOR TXe") 

CALL PLOT4 (2.2HSI) 

WRITE (651200) 

BEtOw PLOTS AEe IMMATURES IN CONTAINERS WITH TxXe 

CALL PLCT1 (32515102111%310) 

CALL PLCT2( GRAPH 6360220 02140020) 

CALL PLOT3 (1H* »sDAY( 1) »SUM2(1),357) 

WRITE (62703) 


FORMAT (1H15*AEs IMMATURES IN CUNTAINERS POSITIVE FOR TXe LARVAE ®* ) 


CALL PLOTS (159 1S5SHAEDES IMMATLRES) 
WRITE (631200) 

BELOW PLOTS TXe LARVAE 

CALL PLOT] (€305159109111,10) 

CALL PLCT2 (GRAPH. 3600200920 020-6) 
CALL PLOT 3( 1H* sDAY(1)25UM3 (1) 2357) 
WRITE (€s704) 

FORMAT (1H1 »*TOXORFYNCEITES LARVAE® ) 
CALL PLCT4 (11ls11HTOX LAR VAE ) 
WRITE (651200) 

FORMAT (1HO253Xe*DAY® ) 

STOP 
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